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ABSTRACT 

The ratio of the Bondi and Jeans lengths is used to develop a cloud-accretion 
model that describes both an inner Bondi-type regime where gas pressure is 
balanced by the gravity of a central star and an outer Jeans-type regime where 
gas pressure is balanced by gas self-gravity. The gas density profile provided 
by this model makes a smooth transition from a wind-type inner solution to a 
Bonnor-Ebert type outer solution. It is shown that high-velocity dust impinging 
on this cloud will tend to pile-up due to having a different velocity profile than 
gas so that the dust-to-gas ratio is substantially enriched above the 1% ISM level. 



1. Introduction 

Laboratory-scale plasma jets have been produced by a magnetohydrodynam ic mecha- 



nism believed analogous t o the mechanism responsible for driving astrophysical jets (iHsu &: Bellan 



20021 : iBellan et al. II2005I ). The laboratory jets are driven by capacitor bank power supplies 
that provide poloidal and toroidal magnetic fields and the jet acceleration mechanism can 
be considered as due to the pressure of the toroidal magnetic field inflating flux surfaces 
associated with the poloidal magnetic field. If these jets are indeed related to astrophysical 
jets, the obvious question arises as to what constitutes the power supply responsible for the 
toroidal and poloidal magnetic fields in the astrophysical situation. Existing models of as- 
trophysical jets typically assume the poloidal field is simply given and that the toroidal field 
results from the rotation of an accretion disk twisting up the assumed poloidal field. The 



- 2 - 



author has developed an ahernate model which postulates that the toroidal and poloidal 
field result instead from a dusty plasma dynamo mechanism that converts the gravitational 
energy of infalling dust grains into an electrical power source that drives poloidal and toroidal 
electric currents creating the respective toroidal and poloidal fields. A brief outline o f how 
infalling charged dust can drive poloidal currents has been presented in iBellanI (120071 ). 



An important requirement of the model is that there should be sufficient infalling dust 
to provide the jet power. It has been well established that the dust-to-gas mass ratio in 
the Interstellar Medium (ISM) is 1%. If one assumes, as has been traditional, that this 
ratio holds throughout the accretion proces s, then the gr avitational energy available from 
infalling dust would be insufficient. However, iFukud (120011 ) has recently shown via numerical 
solution of coupled dust and gas equations of motion that the dust-to-gas ratio can become 
substantially enriched during Bondi-type accretion. 

Star formation has also been previously examined from a molecular cloud physics point 
of view which differs significantly from the Bondi accretion point of view in the treatment of 
self-gravity and inflows. Molecular cloud physics is clearly important because observations 



indicate that stars form within the dense cores of molecular clouds (e.g., see lEvans et al. 



(Ebert 


1955; 


Bonnor 


1956) 



This solution does not take into account flows associated with accretion. Dust dynamics is 
also not taken into account; instead it is typically assumed that the dust-to-gas mass ratio 
is flxed at 1% for all radii. 

To summ arize the above discussion, we note that the Bondi-type analysis reported by 



Fukud (I2OOII ) emphasizes inflow physics and dust-gas couplin g but does no t take into account 
self- gravity whereas molecular cloud analysis such as used by lEvans et al. I ( I2OOII ) emphasizes 
self-gravity, but does not take into account inflow or dust-gas coupling. 

This paper will address the physics of co upled d ust an d gas accretion using a method- 
ology similar in concept to that presented by iFukud (1200 ll ). but extended to bridge the gap 
between Bondi accretion models and molecular cloud force balance models. The enrichment 



mechanism observed by iFukud (1200 ll ) will be examined in detail and will be shown to result 
from the inherently non-uniform nature of the Bondi/molecular cloud/ISM system. Proper- 
ties of dust and gas for radii ranging from the cloud-ISM interface to the inner Bondi region 
will be considered by examining a sequence of successively smaller concentric regions. The 
basic character of each region and its scale, deflned in terms of the nominal radial distance 
from a central object star, is as follows: 



ISM scale: The outermost scale is that of the Interstellar Medium (ISM). The ISM has 
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a gas density ~10^ m~^, a dust-to-gas mass ratio of 1 percent, a gas temperature 
Tg^'^^ ~ 100 K, and is optically thin. The ISM is assumed to be spatially uniform and 
to bound a molecular cloud having radius redge- 

Molecular cloud scale: The molecular cloud scale has much higher density than the ISM 
and is characterized by force balance between gas self-gravity and gas pressure. The 
molecular cloud scale is sub-divided into a large, radially non-uniform low-density outer 
region and a small, approximately uniform, high-density inner core region. Clouds have 
a characteristic scale given by the Jeans length rj. The radial dependence of density 
is provided by the Bonnor-Ebert sphere solution which acts as the outer boundary of 
the Bondi accretion scale. 



Bondi accretion scale: The Bondi accretion scale (jBondilll952l ) is ~ which is sufficiently 
small that gas self-gravity no longer matters so equilibrium is instead obtained by force 
balance between gas pressure and the gravity of a central object assumed to be a star 
having mass M*. The Bondi scale is sub-divided into three concentric radial regions: an 
outermost region where the gas flow is subsonic, a critical transition radius at exactly 
tb where the flow is sonic, and an innermost region where the gas flow is free- falling 
and supersonic. 

Collisionless dusty plasma scale (to he considered in a future 'publication): Free-falling 
dust grains collide with each other in one of the above scales and coagulate to form 
large-radius grains which are collisionless and optically thin. The optically thin dust 
absorbs UV photons from the star, photo-emits electrons and becomes electrically 
charged. The charged dust grains are subject to electromagnetic forces in addition to 
gravity. Motions of charged dust grains relative to electrons result in electric currents 
with a s sociat ed poloidal and toroidal magnetic fields [see preliminary discussion in 



BellanI (120071 )]. 



Jet scale (to he considered in a future publication): The electric currents interact with the 
magnetic fields to produce magnetohydrodyn amic forces whi c h driv e astr ophysical jets 



i n a m anner analogous to that reported in iHsu fc BellanI (120021 ) and iBellan et al. 



The separation-of-scales requirement rj » implies existence of a small parameter 

£bj = — (1) 
rj 

quantifying the separation between the Bondi and Jeans scales. For purposes of relating the 
Bondi and Jeans scales to each other it is convenient to introduce a geometric-mean scale 
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with characteristic length given by 

rgm = -JtbTj = -Je^jTj (2) 

in which case 

Tb < -C rj. (3) 

Gas and dust must be considered separately in each scale. We assume that gas motion 
influences dust motion but not vice-versa so that the dust dynamics can be ascertained after 
gas dynamics has been determined. This assumption is appropriate so long as dust mass and 
energy densities are small compared to corresponding gas densities or if the dust is decoupled 
from the gas. The sequence of scales is summarized in Tables [Hand [21 The numerical values 
of table elements with asterisks will be predicted by the model to be presented here while 
table elements with filled-in numbers represent prescribed physical boundary conditions. 
Here is the gas molecular mass, Ug is the gas density, Tg is the gas temperature, Ug^d 
are the radially inward gas and dust fluid velocities, and pg^d are the gas and dust mass 
densities. Quantities with superscripts 'ISM' are evaluated in the ISM, un-superscripted 
quantities refer to Bondi or molecular cloud regions. 



Our method ology will be conceptually similar to iFukud (l200ll ) but will also differ in 



important ways. iFukud (l200ll ) used a wind equation scheme to model Bondi-type physics 
and, unlike the analysis to be presented here, assumed the Bondi accretion region was directly 
bounded by the ISM, i.e., no molecular cloud region with Jeans- type scaling was taken into 
account. In addition, Fukue assumed that (i) the gas was heated by the combined effects of 
adiabatic compression and friction due to dust-gas collisions, and (ii) the entire region was 
optically thin so that the dust was subject to radiation pressure. Because of the complexity 
introduced by the heating of the gas, Fukue obtained results via numerical solution of three 
coupled differential equations (gas momentum, dust momentum, gas heating). Our approach 
will differ by assuming that (i) the gas is isothermal, (ii) a molecular cloud region lies between 
the Bondi accretion region and the ISM, and (iii) the system is optically thick outside 
the inner part of the Bondi region so that dust is shielded from optical radiation until it 
penetrates to the inner part of the Bondi region. Furthermore, rather than using numerical 
solutions, we will attempt analytic solutions as much as possible, a goal made feasible by 
the isothermal assumption. We believe our assumptions (i)-(iii) reasonably correspond to 
observations that gas is nearly isothermal and that stars form inside the cores of optically 
thick molecular clouds. 
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region 


location 




T 


Ug 


ISM 


r > Tedge 


10^ 


100 K 


* 


B-E sphere 


^gm ^ T <i T^dge 


* 


10 K 


* 


Bondi- Jeans interface 


Tgm = y/rsTj 


* 


)) 


* 


Bondi subsonic 


Tb <r < Tgm 


* 


J) 


< ^KTg/rrig 


Bondi transonic 


r = Tb 


* 


n 


^KTglmg 


Bondi supersonic 


r < Tb 


* 


J) 


> ^jKTglmg 



Table 1: Sequence of regions for gas, * indicates quantity to be discussed/calculated in text 



region 



location 



PdlPg Ud 



ISM 

B-E sphere 

Bondi- Jeans interface r 
Bondi subsonic 
Bondi transonic 
Bondi supersonic 



r > 

y/rBTj <r < 



gm ^ ^JTbTj 
Tb <r < Tgm 

r = Tb 
r < Tb 



0.01 

* 

* 
* 
* 
* 



3 km/s (turbulent acceleration) 
slowing down by gas 
entrained with gas 



Table 2: Sequence of regions for dust physics, * indicates to be discussed/calculated in text 



2. Gas and dust equations 

The steady-state equation of motion for spherically symmetric gas is 



PgUg- 



dUg 

dr 



dPa 



MG 



Pa ^2 +fdrag- 



(4) 



Here Ug is the radial fluid velocity of the gas, Pg = UgnTg is the gas pressure, pg is the gas 
mass density, 

M = M, + [\7ie (P,(0 + PdiO) di (5) 

is the mass inside a sphere of radius r, pd is the dust mass density, and are the mass 
and r adius of the central o b ject s tar, and fdrag is the drag force on gas due to collisions with 
dust. iLamers fc Cassinellil (Il999l ) give the drag force to be 



fdrag = Ud)^^^Crd^ + {Ud - UgY . 

The dust behaves as a zero-pressure fluid, so its steady-state equation of motion is 

dud 



pdUd- 



dr 



MG 

Pd 2 J drag ■ 



(6) 



(7) 



- 6 - 



In the limit th at the integra l in Eq|5]can be ignored, these are the same momentum equations 



considered by iFukud (I2OOII ) 



Observations indicate that the gas temperature of clouds, cores, and accretion disks lies 
in the range 10-30 K indicating that the gas can be considered approximately isothermal. 
If the gas were adiabatic, i.e., if Pg ~ pl^^ then the gas temperature would vary as Tg ~ 
n^/^. Since there is at least a 10'^ increase in gas density from the cloud edge to the Bondi 
region inner free-fall region, an adiabatically compressed gas volume element would have its 
temperature increase at least one-hundred fold as it moved inwards. If in addition, there were 
heating of gas due to dust frictional drag, the gas temperature would increase even more. 
For example, the sound speed in Fukue's Fig. 1 increases by a factor of approximately 40 
from right to left corresponding to a factor of 1600 increase in temperature. Thus Fukue's 
Fig. 1 analysis predicts that the gas temperature increases from 100 K to 1.6 x 10^ K a 
temperature at which the molecular hydrogen would not only have become disassociated 
but would have turned into plasma. Since the observed gas temperature is not 1.6 x 10^ 
K but in fact is clamped at 10-30 K, any heat energy resulting from adiabatic compression 
or friction must be immediately lost, presumably via molecular line emission. In order to 
conform to observations we therefore assume the gas is isothermal, an approximation which 
has the incidental, yet fortuitous side effect of making an analytic approach to the problem 
feasible. 

Since we are assuming that the dust is a perturbation on the gas, we first solve the 
gas equation ignoring the dust, and then use gas equation solutions as inputs for the dust 
equation. This approach is self-consistent so long as the dust is approximately decoupled 
from the gas or if the dust mass and energy densities are much less than the gas mass and 
energy densities. The approach fails in situations where the dust is coupled to the gas and 
attains mass or energy densities comparable to the gas. However, the point at which failure 
occurs is interesting and can be considered a useful prediction of the model. 

We now consider gas dynamics while ignoring the dust in which case the gas equation 
of motion reduces to 



dug 2 d-Pg 



where 




^^=\hr- (9) 



nig 



is the gas thermal velocity. 

Equation [8] has two regimes of interest, namely 
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1. an inner regime where r is so small that ^ An^"^ pg[^)d^ and 

2. an outer regime where r is so large that M^. <^ Ai:^'^ pg{^)d^. 

The former r egime leads t o a Bondi accretion situation characterized by the wind equa- 



tion discussed by iFukud (l200ll ) while the latter leads to a Jeans-scale problem with Bonnor- 
Ebert spheres as the solution. We will examine these regimes separately and then connect 
them using an asymptotic technique. 

Because the configuration is assumed to be steady-state, the gas and dust equations of 
continuity give the respective accretion rates 

Mg = —Anr'^pgUg = const. (10) 

Md = —Aur'^pdUd = const. (11) 
where Ug and Ud are negative, corresponding to radial inward motion. 

3. Gas 

3.1. Bondi accretion region 

In this region it is convenient to normalize the velocity to Cg and lengths to the Bondi 
radius r^, defined as 

= ^ • (12) 

A bar will denote normalized quantities and to avoid confusing minus signs, the normalized 
velocity is defined to be positive so 



Ug 



Ug\/Cg (13) 

f = t/tb- (14) 

Also the gas mass density is normalized to its value at Vgm, the geometric mean of the Bondi 
and Jeans scales (see Eqj2]), i.e., 

Pg = Pg/pgm- (15) 

Thus Pgm is the mass density at a radius much larger than the Bondi length, but much smaller 
than the Jeans length (see Eq|3]). The normalized form of EqElin the Bondi accretion region 
is thus 

=- '^^^Pa - A. (16) 
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The derivative of Eq{TO] gives 

dhipg 1 dug 2 



dr Ug dr r 



(17) 



Combining EqsiT6] and [T7] gives the isothermal wind equation (ILamers fc Cassinellilll999l ) 

i-l - 1) % = (I -I)- (18) 



dr \r r ^ 

Equation [18] gives the condition that Ug = 1 must occur when f = 1 in order for dug/df to 
be non-singular at f = 1. Equation [TSl can be directly integrated to give 

V? 2 

-j- -\YiUg = 2\Tif + - + A (19) 

where A is a constant to be determined. The sonic condition = 1 at f = 1 gives A = —3/2 
so EqJT9] becomes 

ugf^ - exp (I - ^ + = 0. (20) 

Equation [201 the solution to Eq{THl is a transcendental expression relating Ug and f and 
having two distinct roots. One root has Ug a monotonically increasing function of f and one 
root has Ug a monotonically decreasing function. Both roots have -Ug = 1 at f = 1 . The 
root where Ug is a monotonically increasing function of f (subsonic at small f, supersonic 
at large f) is relevant to stellar winds, while the monotonically decreasing root (supersonic 
at small f, subsonic at large f) is relevant to the accretion problem discussed here and by 
Fukue. Since Ug is subsonic at large f and supersonic at small f, Eq{T9] has the following 
limiting forms 



luMg = ln2 Inr for r^l (21) 



3 

Inw,, = - 21nf + - for f > 1. (22) 

The change in the slope of In-u^ as a function of Inf from —1/2 at small f to —2 at large 
f is evident in Fig. 1 of Fukue. 

The mass accretion rate can be determined by evaluation at f ^ 1 since in this limit 
Eqj20] gives 

-ex (-- - 



""s ^ :^*^^p(~~+o) r > 1. (23) 



Using Eq|23] in Eq{TO] and then evaluating at Tgm gives the Bondi mass accretion rate 

Mg = ,!^^W^ (24) 
r 

^9 
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where 



r]{eBj) = exp ( -2^^ + - ) (25) 



is a dimensionless quantity of order unity which occurs repeatedly in the analysis to follow. 
The radial dependence of the gas mass density in the Bondi regime is given by 

Pg becomes independent of r in the large r limit of the Bondi region since in this limit Eql23] 
shows that Ug ~ This independence permits matching the large r limit of the Bondi 

region solution to the small r limit of the molecular cloud region solution (given in Sec J3.2l 
below) since both these solutions are independent of r. 



In the small r limit, Eq. [19] gives 



u 



9 



2/fi/2 ^27) 



which corresponds to the free-fall velocity. Combining this with EqsiTO] and [24] gives the 
mass density for r << to be 

= (28) 

The enclosed gas mass at is 

r-TB /■! 4 

/ A-n-r'^pgdr ~ 47rpg„r| / f'^pgdf = -Tiripgmr% . (29) 
JQ Jo 3 



3.2. Molecular cloud region 

Again ignoring dust mass and dust drag, but now assuming f is sufficiently large that 
M, < J"^ An^^ (pg(0) d^, EqE] becomes 

+ ^ r 4^^2 (^)^^ ^0 (30) 
Pg dr Jj.^ 

where the subsonic condition << has been used. Gas pressure is now balanced by 
self-gravity rather than by the central object gravity. Equation [50] can be recast in standard 
form as 

cl 1 d f r"^ dpA , ^ 
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We now define a new dimensionless lengtli suitable for Eqj3T| namely 



r 

f = — (32) 
rj 



where 

r T = 



= -7==^ (33) 



is the Jeans length. 

Equation [31] assumes the dimensionless form 

1 d f r"^ dpg 



f 2 dr \ pg df 



+ -pg =0. (34) 



For small f, we wish to have a solution which is approximately constant in order to be 
consistent with observations. This constant should be unity in order to connect with the outer 
limit of the Bondi solution since the outer limit of the Bondi solution has been constructed 
so as to give p pgm for r = rgm- An approximate solution to Eql3l] that does this is 
Pg = l/(l+f^/6) but this solution fails at large f. However, Equation[3l]has the exact solution 
Pg = which is valid at all f. Thus, a solution which behaves as Pg = 1/(1 + f^/6) for 
small f and as = for large r would do the job provided these sub-solutions smoothly 
merge into each other and do not introduce any spurious effects at intermediate values of r. 
The Fade approximation 

Pgir) = ^ ^ ox (35) 

agrees with a direct numerical solution of EqlM] with an error less than 20% in the range 
< f < 9, an accuracy which is more than adequate for the present analysis. The factor 
(f^/(15 + r^))^ can be thought of as a 'switch' which causes to become f^/2 at large f^. 

At Tgm, the value of f is f = rgm/rj = ^Ebj and so EqJMl matches to the Bondi solution 
at Tgra since by assumption ebj ^ 1- 

The combined solution spanning the range < r < 9rj is therefore 

for r < r, 



4r2|Mg(r)|c3 Pgm 



1+2 



6r T \ \ 15r j-\~r 



for rqm<r < 9r j 



gm 

(36) 



where Ug{r) is given by Eq|20]and both outer and inner solutions converge to pgm at r = rgm- 
Figure [T]( a) plots %(r) = pg{r)/mg using Eql36] with the assumptions that Sbj = 0.3, 
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M, = OAMq, Tg = 10 K, T/^^^ = 100 K, and ^ = lO^m'^. The values of r^, Vgm, and 
rj are indicated by vertical dashed lines. The solid line in Fig. [^b) plots the corresponding 
gas velocity Ug{r)] note that Ug = Cg at rs and that the gas is subsonic to the right of r^. 

Specification of Ebj determines pgm since substitution of Eqs{T2] and [33] into Eq{T] gives 

Pa^ = 'BJ^^^^- (37) 

The separation-of-scales condition is consistent with the assumption -C /J^ 4:'7i^'^pg{^)d^ 
for r>rj since /J"/ "^^^^PgiOd^ - ^T^pgrnr]/^ = MJQebj. 

Equation [371 can be used with Eql29]to calculate the total gas mass enclosed at to 

be 



A7rr^Pgdr= |4j^* (38) 
which is consistent with the Bondi region assumption that f^^ Airr'^pgdr <^ M^. 

Equations [21] and [37| show that the gas mass accretion rate can be expressed as 

Mg= |4j^. (39) 



The gas fluid velocity at rgm is 



M, 



Ug (rgm) = — — = -EBjVCg. (40) 



am 



Bondi accretion will presumably increase M*. Because of the existence of jets, not all 
accreting matter will do this. However, for purposes of estimation, if one assumes that all 
accreting material causes an increase in M^,, then tb which is proportional to will increase. 
On the other hand rj does not depend on M,, and so will remain constant. Thus, tb/tj will 
be a slowly increasing function of time and eventually the presumption that tb '^tj fails. 



3.3. Interface between molecular cloud region and the ISM 

Although gas density and temperature in reality change gradually on entering the cloud 
from the ISM, for simplicity we assume here that these changes occur in a narrow layer. 
Integration of EqH] across this layer while taking into account mass conservation gives con- 
tinuity of PgU^g + Pg across the layer and integration of the mass conservation equation gives 



3 

cr 
3 



10" 
10000 
in 1000 
100 
10 

1 

10.00000 
1 .00000 



0.00100 

0.00010 
0.00001 

1 .0000 
0.1000 

5 0.0100 
I. 

0.0010 
0.0001 



10' 



10^ 



gm 



gm 



gm 



10'^ r. 



gm 

m 

10" 
ou 



(d) 



10' 



' edge 









(b) I 








/ p 





edge 









(c) \ 











edge 



edge 



10=" 



Fig. 1.— Model predictions assuming Ebj = 0.3, M, = OAMq, T,, = 10 K, T/^^ = 100 K, 
and Ug^^^ = 10^m~^. (a) gas density, (b) solid line gives gas velocity, dashed line gives dust 
velocity, (c) dust density, (d) dust-to-gas mass ratio. 
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continuity of pgUg across the layer. This yields Rankine-Hugoniot equations of the form 

n,(re,,e) (1 + {ug{r^.,e)f) = nf ( Tf ^V^, + (^^0') (41) 

%(re.,e)%(re.,e) = ^ (42) 

where ng^'^,TisM, and Tg are specified and u^^'^^ ,redge, Ug^r^dge), and ng{redge) are to be 
determined. The bar indicates that densities are normalized to and velocities are nor- 
malized to Cg = ^y nTg/nig. Using Eq|42]to substitute for ng{redge)ug{redge) in EqjH] gives 

1 1 \ n^iSM 



However using EqJTO] evaluated in the ISM at rg^^g, Eql39]and Eqj33] give 



r 



2 



" = -Velj ^ (44) 

while Eqj35] gives 



edge 



^^(redge) ^ (45) 
edge 



SO eliminating r"^]/ ridge between these equations gives 



^ISM-ISM ^ _!^^^(,^,^^). (46) 



2 

Substituting for nl^^^ui^'^^ in EqHSl gives 



9 9 



a quadratic equation in ngir^dge)- Since e^j is assumed small, an approximate solution to 
EqJTTl may be obtained by balancing the last two terms in which case 

^/SM rpISM 



Equation US] then gives 



so using Eqj33] 



J - 9 9 ^^Of) 



idge (1 + Tg 



2 , 1+ 4jr/74 
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Equation HH gives 

rpISM 



..ISM _ 2 _g "-g /c-.N 

~ ^^^^ 2T, ii + e%jvy^y ^ ' 



which combined with Eqj42] and Eqj48] give 



%(re.,e) = - . (52) 



The gas mass between r j and Vedge is 



edge 



1 + e%jvy^ 



The mass inside the core region of radius r j is using EqsJ35l 1331 and [37] 

Table [3] shows quantities predicted by this analysis for Sbj = 0.3, = O.4M0, T^, = 10 
K, Tg^^^ = 100 K, molecular hydrogen, and Ug^'^^ = 10^ m~^. Note that the gas flow velocities 
are subsonic (i.e. much smaller than the random or thermal velocity) both in the portion of 
the cloud external to and in the ISM since in the portion of the cloud external to gas 
flow velocities Ug are much slower than Cg = ^JKTg/mg = 200 m/s and in the ISM the gas 
flow velocity u^'^*^ is much slower than c^"^*^ = ^ kT^^^^ /nig = 640 m/s. 



4. Dust 



4.1. Overview of dust 



(LiUey 


1955: 


Zubko et al. 


2004) 



1% ratio also holds in a dense molecular cloud. However, iPadoan et al. I (120061 ) have recently 
cast doubt on such an assumption arguing that if the dust-to-gas mass ratio were indeed 
1% in a dense molecular cloud, then the dust and column densities should have the same 
spectral power law, i.e., the spatial Fourier power spectra S{k) o f dust and gas i r nages should 



be characterized by the same power law S{k) ~ In fact, iPadoan et al. I (120061 ) found 

that p for the dust differs from p for the gas, indicating a lack of local proportionality between 
dust and gas. 
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1 

value 




units 


1 

value 


units 


Keierence 




6.5 X 


10^^ 


m 


4.3 X 10^ 


au 


EqiH 




1.2 X 


10^^ 


m 


7.9 X 10^ 


au 


Eql2] 


rj 


2.2 X 


10^^ 


m 


1.4 X 10"^ 


au 


Eql33] 




3.1 X 


10^ 




3.1 X 10^ 


cm~3 


EqET] 




-91 




m/s 






Eqgn] 


Cg 


2.0 X 


10^ 


m/s 






EqE] 


Mg 


1.7 X 


10^^ 


kg /s 


2.7 X 10"' 


Mo yr-i 


EqlSi 


"^edge 


1.7 X 


10^^ 


m 


1.1 X 10^ 


au 


Eql50] 


^edge / ^ Jeans 


( . ( 












^g ('"edge) 


-14 




m/s 






Eqisa 




-1.4 


X 10^ 


m/s 






EqED 


„ISM 


6.4 X 


102 


m/s 






Eqia 




2.1 X 


1031 


kg 


10 




Eql53] 




4.4 X 


1029 


kg 


0.2 




EqEl 


^r<rB 


1.8 X 


1028 


kg 


0.01 


Mo 


Eql^l] 



Table 3: Calculated quantities for e^jj =0.3, M, = O.4M0, Tg=10 K, = 2 x 10^ m/s 



Takeuchi et al. I (120051 ) have also pointed out that the dust/gas mass ratio in a cloud 



need n ot be the same as in the ISM because gas and dust might be subject to different pro- 



cesses. 



Goldsmith et al. I (119971 ) showed in their Fig. 22 that the observed dust-to-gas ratio 
varies by an order of magnitude as a function of position in a system of molecular clouds, i.e., 
the observed dust-to-gas mass r atio in a mole cular cloud system is not, as in the ISM, fixed 
at 1%. As mentioned earlier, iFukud (I2OOII ) presented numerical calculations specifically 
demon s tratin g that accretion can enrich the dust-to-gas ratio. As also mentioned earlier, 
Fukud (I2OOII ) assumed that a Bondi accretion region was directly bounded by the ISM, i.e., 
did not take into account the molecular cloud scale where gas self-gravity is important. This 
omission resulted in calculated mass accretion rates lower than typical observed values (e.g., 
Fukue's nominal accretion rate was Mg ~ lO^^^ Mq yi~^ compared to the observed nominal 
M„ ~ 10~9 Mq to 10~^ Mq rates for pre-main sequence YSO's reported by lHartmann et al. 



(119981 )). Since gas and dust density scale linearly with accretion rate, the low nominal accre- 
tion rate in Fukue's calculation corresponded to a very low dust density and hence optically 
thin dust. Higher accretion rates give higher dust densities and optically thick dust densities 
that correspond to the observed opacities of clouds. Nevertheless, Fukue's calculation had 
the interesting feature of predicting that in certain circumstances the accretion process could 
cause a thirty-fold increase in the dust-to-gas mass ratio. 
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We now argue that this enrichment reported by Fukue should be a basic property of 
dust /gas accretion systems, so the dust/gas mass density ratio in a cloud core should in 
general be substantially enriched relative to its 1% ISM ratio. Because the accretion region 
is now assumed to be bounded by a cloud core rather than by the ISM, this enrichment effect 
is demonstrated here in association with more realistic mass accretion rates. 

Collision properties of gas and dust vary considerably in ISM, cloud, cloud core, and 
accretion regions. This variation is because the scale lengths are different and because the 
mean free paths are different. Dust is coUisionally decoupled from gas in some regions, but 
strongly coupled in others. In particular, it turns out that dust is decoupled from gas in the 
ISM, strongly coupled in the molecular cloud and cloud core, but then can become decoupled 
again when dust-dust collisions result in dust coagulation. Dust-dust coUisionality is closely 
related to dust opacity since the mean free path for dust-dust collisions is 1^ = ^/nd<7d while 
the distance for unity optical depth is the same if light scattering is geometric, but larger 
by the factor 1/Qext, where the efficiency factor Qext accounts for weaker than geometric 
scattering (i.e., Rayleigh) and for absorption. 

Let us start with the ISM. Because of the low gas density in the ISM, dust is coUisionally 
decoupled from gas so if the dust becomes charged, electromagnetic mechanisms can acceler- 
ate ISM dust grains t o velocities much higher th an the ~ 640 r a/s IS M gas thermal velocity. 



Mever et al. I fll998f ). lYan fc LazarianI ([20031), ISlavin et al. I mm . and IShukla fc Stenflo 



(120051 ) have provided examples of such dust acceleration mechanisms. 



In contrast, dust is coUisionally coupled to gas in molecular clouds because of the high 
gas density. In cloud cores and/or accretion regions the dust density becomes so high that 
dust-dust collisions also become important. The condition for dust-dust collisions to be 
important corresponds approximately to the condition that the dust is optically thick. Thus, 
if there is an inner region where dust-dust collisions are important, regions external to this 
inner region will be shielded from central object optical emission and so will not experience 
any radiation pressure from the central object. Dust-dust collisions result in dust coagulation, 
a condition where the dust radius increases while the overall dust mass density remains 
invariant. Coagulation has the dual effect of causing the dust to revert to being optically thin 
and collisionless. Collisionless, coagulated dust exposed to stellar radiation becomes charged 
via photo-emission and so must be described by dusty plasma dynamics rather than by gas 
dynamics. This dusty p lasma regime will be considered in a future publication; preliminary 
results are presented in iBellanI (120071 ). 
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4.2. Dust grain size assumption 



Mathis, Rumpl, and Nordsieck (1977) reported observations showing tliat ISM dust 
grains liave a size distribution scaling as r^^'^ witli a nominal lower limit radius riow = 0.005 
/im and a nominal upper limit Vhigh = 0.25 /im. This size dis t ribut ion is commonly referred 
to as the MRN dust size distribution. iMiyake &: Nakagawa] (Il993| ) proposed that the dust 
size distribution in a circumstellar disk will have larger grains than given by the MRN 
distribution because in a circumstellar disk where the dust density is hig h, dust coagulates 



due to dust-dust collisions. More recently, iPrzygodda et al. I ( l2003l ) and Ivan Boekel et al. 



( I2OO3I ) have report ed direct observational eviden ce of grain growth in circumstellar disks 
while, in addition, iDuUemond fc DominikI (120051 ) provided detailed calculations showing a 
strong tendency for dust grain growth when dust grains collide with each other. 

The r"^-^ dust size distribution has important well-known statistical properties. Specifi- 
cally, an r^^'^ distribution implies (i) most of the mass is contained in the very small number 
of heav y grains and (ii) most of the surface area is contained in the very large number of light 
grains (ISmith fc Dweklll998l ). Although the lightest dust grains provide most of the surface 
area, their radii are so much smaller than visible light wavelengths that their scattering is 
Rayleigh rather than geometric. Rayleigh scattering is weaker than geometric scattering by 
a factor r^/ where A is the wavelength of the radiation. Hence, despite being numerous, 
the very small radius grains are ineffectual at contributing to the optical depth at visible 
wavelengths so visible scattering will be due mainly to the larger grains which also contain 
most of the mass. Hence, for purposes of both mass inventory and optical depth we will 
consider that dust grains have a nominal radius = 0.1 /im. 



4.3. Dust velocity in the ISM 



If the velocity of individual dust grains in the ISM were the result of thermal equilibra- 
tion with gas molecules, the dust kinetic temperature would be the same as the gas molecule 
kinetic temperature in which case the random velocity of dust grains in the ISM would be 
y/5Af ^ ^ KTg/nid = 0.01 m/s. This velocity is so small that there would be insufficient 
dust flux entering a molecular cloud to populate the cloud with dust in a reasonable time. 
It has been proposed by several authors that dust attains a much larger random velocity 
in the ISM because dust is nearly coUisionless in the ISM so that collective coUisionless 
mechanisms such as shocks or turbulen c e could accelerat e dust t o very large y e locitie s. As 



mentioned earlier , iMeyer et al. I (1l998l ). lYan &: Lazariaru ((20031), ISlavin et al. I (1200J), and 



Shukla &: Stenflol (120051 ) have presented coUisionless mechanisms whereby dust grains are 
accelerated to velocities in the 10 km/s range. 
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A typical argument for coUisionless acceleration is that a spatially and temporally ran- 
dom magnetic field develops with energy density in thermodynamic equipartition with the 
gas, i.e., 5V2/io ~ ^^kT/^^'^ so 



B ~ *^kT/^*^ ~ 2 X 10"^° T. (55) 

Gradients in this microgauss magnetic field then accelerate charged dust grains via electro- 
magnetic forces so that the dust random kinetic energy becomes comparable to the magnetic 
energy density, i.e., 

U'''{Ur'T-f- (56) 

where U^^'^ is a random velocity that could be pointing in any direction. The nominal 
random velocity f/j"^*^ attained by the accelerated grains is thus the dust Alfven velocity 
given by v\ = B"^ / ^^p^^^^^ . Combining Eqsl55] and 1561 gives the random velocity to be 



V ^9 Pd 

a relationship dependent on the dust-to-gas mass density ratio. Equation EH gives U^^^^ ~ 
9.3 km/s for a nominal 1% ISM dust-to-gas mass ratio and T^^^^ = 100 K. This velocity 
is six orders of magnitude larger than the velocity predicted by thermal equilibration of 
individual dust grains with gas molecules. Because the dust is assumed to be accelerated by 
coUisionless processes, it will not in general have a Maxwellian distribution, and in fact it 
would be reasonable to expect that each dust grain has a velocity ~ U^^^'^ but with a random 
direction. This situation would be similar to the neutrons emitted from a fission reaction 
since such neutrons all have the same kinetic energy but have random velocity directions. 

The mean radial velocity of dust grains entering a spherical molecular cloud can be 
determined by considering the radial velocity component of these dust grains at the surface 
of the sphere. If s denotes the inwards normal to the surface at some point on the surface 
of the sphere, then the radial velocity of a dust grain entering the sphere at this point is 
Uji^^ cos 6 where 6 is the angle between the dust grain vector velocity and s. The angle 
9 must lie between — 7r/2 and 7r/2 since dust grains having 9 outside this range will not 
enter the sphere. Other than being restricted to the range — 7r/2 < 9 < 7r/2, the angle 9 is 
random. The dust grains have equal probability of having any angle 9 between and 2tt so 
the probability of a dust grain angle being between 9 and ^ -|- c?^ is d9/2TT. The mean radial 
velocity of dust grains entering a molecular cloud from the ISM is thus 

u'SM ^ r" ^isM ^ — = -U'/^' ^ 3 km/s. (58) 

J-7T/2 271 TV 

The question of whether dust exits a molecular cloud will be addressed in the next section. 
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4.4. Dust entrainment by gas in molecular cloud 

The dust equation of motion can be expressed as 



du G ( 



(59) 

where it is recalled that both Ud and Ug are negative corresponding to radial inward motion. 

We wish to investigate the possibility that Pd/ Pg could increase from its 1% ISM value 
as dust and gas enter the cloud, core, and Bondi accretion region. However, we will assume 
that Pd/ Pg remains sufficiently small compared to unity that the back- reaction of dust on 
gas dynamics can be ignored. Thus, Ug{r) and Pg{r) are assumed to remain as discussed in 
Sec J3.1l and 13. 2[ Dust enters the cloud from the ISM with a velocity much higher than the 
gas thermal and fluid velocities but then collides with the gas, thereby slowing down until 
moving at nearly the same inward flui d velocity as the gas. Collisions thus cause the dust 



to become entrained by the gas inflow (lFukuell200ll ) 



The coUisional slowing down can be estimated by assuming that the gravitational force 
term in Eql59] can be ignored because of the high dust velocity. Furthermore, Ud is much 
larger than both Ug and Cg when the dust enters the cloud from the ISM. In this case and 
using 

o-d 3 



rrid ^Tdp'd* 

where pj/^* ~ 2 x 10"^ kg m~^ is the intrinsic mass density of the dust, EqlSIJ] becomes 



(60) 



dud _ _ 3pg(r 
dr '^TdP^d'^ ' 



T—M- (61) 



Integration of Eq. [6T] starting from the cloud edge and going inwards gives 

Udir) = Udir^dge) exp (- i^^^) " (^2) 

The dashed line in Fig. W(Jo) shows the dust velocity given by Eql62]for the gas density given 
in Fig. [H^a) and assuming a 3 km/s dust mean radial entrance velocity. It is seen that the 
dust rapidly slows down until it achieves the same velocity as the gas at which point it is 
assumed to be entrained by the gas. Dust will thus become entrained by gas at a radius 

"^entrain whcrC Tgntram SatisflcS 



'edge 3pg(^r) 

Ardpf 



-f"-p(-/ T^tdr -uf^^r.Cg/20. (63) 
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Using Eq|58]to estimate u^/'^ , Eqj63] becomes 

Mi) . ,„ f I ^ 5.7 . (64) 



If the gas density in the r » rj portion of the cloud is sufficiently large to provide en- 
trainment, then Pg{r) can be approximated in this region as Pg{r-) ~ 2pgmr']/r'^ so Eqj6^ 
becomes 

1^ ( ^ V 5.7 (65) 

^^dPd V '^entrain ^edge / 

Table 3 gives r^dge/rj ~ 7.7 and since ^Pg^rj /2p''^^rd = 1.7 x 10^ this gives r entrain - 6rj. 
As assumed, the dust becomes entrained by the gas at a location in the r » rj portion of 
the cloud. 

The entrainment of dust by gas means that the radial flow of dust from the ISM to the 
cloud is one-way inwards, i.e., there is no radial outward flow of individual dust grains from 
the cloud back to the ISM. This one-way inward behavior of dust grains is in contrast to 
gas. Because gas molecules collide with each other, some gas molecules have inwards radial 
velocities and some have outwards radial velocities. However, there are more inwards than 
outwards moving gas molecules and the net inwards gas fluid velocity at r^dge is a consequence 
of this difference. The situation is somewhat analogous to a home vacuum cleaner ingesting 
macroscopic particles and air; individual air molecules in the vacuum hose have inwards or 
outwards velocities with magnitude of the order of the air thermal velocity, the mean air 
velocity in the hose is inwards and much slower than the air thermal velocity, and all the 
macroscopic particles in the hose have inwards velocities. 



4.5. Enrichment of dust-to-gas mass ratio 
Equations [10] and [11] can be evaluated at any radius. Evaluating at r^dge and at 



gm 



gives 



4vrr,2,,,pf *^«f = ^7irl^Pg{rgn.)u^irgm) (66) 
4vrr,%,pf = ^7Trl^Pd{rgn,)u,{rgm). (67) 



Dividing EqES]by EqEZ] gives 



gm ) 



Since the dust is certainly fully entrained by the time it reaches we may assume 
Udi^gm) /ug{rgm) = 1. Using Eq(5T] to give m^'^*^ and Eqs. [57] and [58] to give u^^^'^^ Eql68] 
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becomes 

PAr^m) ^ 2 + e%jr,^l2 / 2T, 

Fig. [Uj^c) shows the dust density and Fig{T]^d) shows the dust-to-gas mass ratio as function 
of radius, showing that enrichment occurs so that the dust-to-gas ratio increases from its 
0.01 ISM value to become somewhat less than unity. For the example parameters in Table 3 
and a 1% ISM dust-to-gas mass ratio, Eqj69] predicts PdiTgm) / PgiTgm) = 0.2 for Ebj = 0.3, 
Tg = 10 K, and T^^'^^ = 100 K; this corresponds to the saturated value in Fig{T](d) and 
represents a twenty-fold enrichment of the dust-to-gas mass ratio over its ISM value. 

two-lane highway merging one-lane highway 




truck speed twice car speed truck speed same as car speed 

truck density half car density truck density enriched 2x 



Fig. 2. — Because trucks entering the one-lane segment slow down whereas cars do not, the 
truck-to-car density ratio is enriched in the one-lane segment compared to the density ratio 
in the two- lane segment. 



An everyday example of the enrichment effect predicted by Eql68] is where traffic on a 
two-lane highway merges to flow on a one-line highway as sketched in FigO In the two-lane 
segment (left of figure), the trucks move at twice the speed of the cars and have half the 
density of the cars. The cars are assumed to have the same speed in the two-lane and one- 
lane highway segments, i.e., the cars do not slow down. The trucks have to slow down when 
they enter the one-lane segment because in the one-lane segment the trucks must go at the 
same speed as the cars. Slowing down by a factor of two causes the number of trucks per 
length of road to double because in steady-state the fiux of trucks (product of the truck speed 
and the trucks per length of road) must be the same in the one- and two-lane segments. The 
density of trucks to cars is thus enriched by a factor of two in the one-lane highway segment 
compared to its value in the two-lane highway segment. 
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4.6. Radiation pressure and optical depth 



Abbas et al. I (120031 ) used a calibrated electrodynamic balance to measure the radiation 
pressure force exerted on a =0.125 /im particle by a 0.532 nm wavelength laser. They 
found that the force exerted by a laser intensity 8 x 10^ W was 3 x 10"^'^ N. If such a 



particle were in the gravitational field of a star with luminosity 
force on the particle would thus be 

Ue'-^'-) 3 X 10-17 

i'rad = — : z rr:^ iN 



47rr^ 



X 105 



the radiation pressure 



(70) 



where 



r r 



Qextnd{C)(^ddi 



(71) 



is the optical depth for light traveling from the star to radius r and Qext is the extinction 
efficiency. The gravitational force would be 



(72) 



SO the ratio would be 



rad 



L,e-^W 3 X 10-17 1 
47r 8 X 105 m„M,.G' 



Assuming nig = 10 kg for a nominal = 10 m dust grain gives 



rad 



0.9e-"(") 







(73) 



(74) 



so a dust grain could in princi ple be subject to sig nificant radiation force if r <C 1 as was 
proposed by Fukue. Table 1 of lAbbas et al. I ([20031) gives the calculated radiation pressure 



efficiency Qpr = 0.28 for a, = 0.125 /xm dust grain illuminate d by 532 nm l ight. This 



is in reasonable agreement with their measured value Qpr = 0.22; lAbbas et al. I (120031 ) also 
presented a calculated extinction efficiency Qext = 0.33 and a calculated scattering efficiency 
Qsca = 0.16 that together would be in reasonable agreement with the measured Qpr. Figure 
[U^c) shows that dust for the representative parameters discussed here has a density ~ 1 
m~^ at tb and that this density increases for r < r^. Table 3 gives = 6.5 x 10^^ m. Thus, 
a lower bound for the optical depth at can be estimated using 



r{rB) > Qextnd{rB)(TdrB = 0.33 x 1 x tt x (IQ-^ x 6.5 x 10^^ = 6.7 



(75) 



Star light would thus be attenuated by a factor greater than exp(— 6.7) ~ 10 and so optical 
radiation from the star should be fully extinguished by the time it reaches in which case 
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there would be no significant radiation pressure on dust at radii of the order of or larger. 
The reasons why we differ from Fukue by concluding that radiation pressure is unimportant 
are (i) Fukue's adiabatic assumption resulted in gas that was 1600 times hotter at small r so 
for a given pressure the gas density would be 1600 times less and so for a given dust-to-gas 
ratio, the dust density wou ld be much lo A ver, ( ii) Fukue also used a much lower value of 
Qext than the value given in lAbbas et al. I (120031 ) . 



5. Discussion 

The logical requirement that the Bondi length must be smaller than the Jeans length 
introduces the dimensionless parameter Ebj given by the ratio of these lengths. An inter- 
mediate distance Tgm can be defined as the geometric mean of the Bondi and Jeans lengths 
so, if the two lengths are well separated, Tgm is much larger than the Bondi length but much 
smaller than the Jeans length. Thus, r^^ constitutes the large r limit of the solution to the 
Bondi problem and the small r limit of the solution to the Jeans problem. Matching the 
Bondi and Jeans solutions at Tgm provides a solution valid over scales ranging from smaller 
than the Bondi length to larger than the Jeans length. This allows accounting for transonic 
fiow, mass accretion, and gravitational self- confinement in one self-consistent solution to the 
gas equation. 

By assuming that the dust mass density and energy density are not larger than the 
corresponding gas densities, the gas equations can be evaluated ignoring interaction with 
dust. Once the gas behavior has been worked out, the dust can be considered as moving 
through a pre-determined gas profile that retards the dust due to frictional drag. The 
dust slows down greatly whereas the gas velocity changes relatively little, a distinction that 
enriches the dust-to-gas mass ratio. Plausible values of Ebj suggest that the dust-to-gas 
mass ratio will always be enriched to be somewhat less than unity, because at unity the 
assumption that the gas is unaffected by dust fails. The high dust density at small radius 
means that dust-dust collisions will become important and cause coagulation of dust. This 
coagulation will cause the dust to become optically thin and colhsionless again. Because the 
coagulated dust is optically thin, it will absorb photons from the star and become electrically 
charged due to photo-emission of electrons. Two-dimensional motion of colhsionless, c harged 



dust i n a gravitational field resulting in a poloidal field dynamo has been discussed in iBellan 



(120071 ): three dimensional motion and resulting poloidal/toroidal field dynamo action will be 



discussed in a future publication. 



The author wishes to thank an anonymous referee for his/her thoughtful comments. 
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